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Abstract. An improved version of the Cascade-Exciton Model (CEM) of nuclear reactions 
realized in the code CEM2k and the Los Alamos version of the Quark-Gluon String Model 
(LAQGSM) have been developed recently at LANL to describe reactions induced by particles 
and nuclei for a number of applications. Our CEM2k and LAQGSM merged with the GEM2 
^ ' evaporation/fission code by Furihata have predictive powers comparable to other modern codes 

\ and describe many reactions better than other codes; therefore both our codes can be used as 

■ reliable event generators in transport codes for applications. During the last year, we have made 

' a significant improvements to the intranuclear cascade parts of CEM2k and LAQGSM, and have 

'n— «' , extended LAQGSM to describe photonuclear reactions at energies to 10 GeV and higher. We 

have produced in this way improved versions of our codes, CEM03.01 and LAQGSM03.01. 
For special studies, we have also merged our two codes with the GEMINI code by Charity 
and with the SMM code of Botvina. We present a brief description of our codes and show 
illustrative results obtained with CEM03.01 and LAQGSM03.01 for different reactions compared 
with predictions by other models, as well as examples of using our codes as modeling tools for 
nuclear applications. 
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^ . 1. Introduction 

^ \ Following an increased interest in nuclear data for such projects as the Accelerator 

Transmutation of nuclear Wastes ( ATW) , Accelerator Production of Tritium (APT) , Spallation 
Neutron Source (SNS), Rare Isotope Accelerator (RIA), Proton Radiography (PRAD) as a 
radiographic probe for the Advanced Hydro-test Facility and others, for several years the US 
Department of Energy has supported our work on the development of improved versions of the 
Cascade-Exciton Model (CEM) of nuclear reactions j^EllH] and the Los Alamos version of the 
Quark-Gluon String Model (LAQGSM) ^ to describe reactions induced by particles and nuclei 
at energies up to about 1 TeV/nucleon. To describe fission and production of light fragments 
heavier than ^He, we have merged both our codes with several evaporation/fission/fragmentation 
models, including the Generalized Evaporation/fission Model code GEM2 by Furihata j^. Our 
codes perform as well as and often better than other current models in describing a large variety 
of spallation, fission, and fragmentation reactions, therefore they are used as event-generators 
in several transport codes. The status of our codes as of the middle of 2004 together with 
illustrative results and comparisons with other models can be found in [SJ [7| [H] and references 
therein. Here, we present several improvements developed during the last year leading to the 
new versions of our codes, CEM03.01 and LAQGSM03.01 01^]. 



2. Basic assumptions of CEM and LAQGSM 

The Cascade- Exciton Model (CEM) of nuclear reactions was proposed initially at the Laboratory 
of Theoretical Physics, JINR, Dubna JOl to describe intermediate-energy spallation reactions 
induced by nucleons and pious. It is based on the Dubna IntraNuclear Cascade (INC) [H^ and 
the Modified Exciton Model (MEM) |12| I13j . It was extended later to consider photonuclear 
reactions and to describe fission cross sections using different options for nuclear masses, 
fission barriers, and level densities, and its 1995 version CEM95 was released to the public 
via NEA/OECD, Paris and RSICC, Oak Ridge (see corresponding references in |3]). 

LAQGSM is an extension by Gudima, Mashnik and Sierk of the Quark Gluon-String Model 
(QGSM) by Amelin, Gudima and Toneev done at Los Alamos |4j and is intended to describe 
both particle- and nucleus- induced reactions at energies up to about 1 TeV/nucleon. LAQGSM 
is based on a time-dependent version of the intranuclear-cascade model (different from the one 
used in CEM) developed at JINR, Dubna, often referred to in the literature simply as the Dubna 
intranuclear Cascade Model (DCM) (see [TSl and references therein). LAQGSM [311^] differs 
from QGSM p3j by using an extended and significantly improved version of the INC model 
mini) by replacing the preequilibrium and evaporation parts of QGSM described according to 
the standard CEM ^U] with the new physics from CEM03.01 and also has a number of 
improvements and refinements in the Fermi break-up and coalescence models in comparison 
with QGSM jl4j. 

CEM and LAQGSM assume that the reactions occur generally in three stages. The first 
stage is the INC, in which primary particles can be re-scattered and produce secondary particles 
several times prior to absorption by, or escape from the nucleus. When the cascade stage of 
a reaction is completed, CEM and LAQGSM use the coalescence model as described in JH] 
to "create" high-energy d, t, ^He, and ^He by final-state interactions among emitted cascade 
nucleons. The emission of the cascade particles determines the particle- hole configuration, Z, A, 
and the excitation energy that is the starting point for the second, preequilibrium stage of the 
reaction. The subsequent relaxation of the nuclear excitation is treated in terms of an improved 
version of the modified exciton model of preequilibrium decay followed by the equilibrium 
evaporation/fission stage of the reaction. Generally, all four components may contribute to 
experimentally measured particle spectra and other distributions. But if the residual nuclei 
after the INC have atomic numbers with A < 12, CEM03.01 and LAQGSM03.01 use the Fermi 
break-up model JB] to calculate their further disintegration instead of using the preequilibrium 
and evaporation models. Fermi break-up is much faster to calculate and gives results very similar 
to the continuation of the more detailed models to much lighter nuclei. 



3. Recent developments in CEM03.01 and LAQGSM03.01 and illustrative results 

First, during 2004 we incorporated into LAQGSM the improved approximations for the total 
elastic and inelastic cross sections of hadron-hadron and photo-hadron elementary interactions 
developed previously for the code CEM97 dH and CEM2k P (see details in |T7j^. 

Second, the double differential distributions of secondary particles from elementary NN and 
jN interactions were simulated by CEM2k (and all its precursors, as well as by LAQGSM and 
its precursors at energies below 4.5 GeV/A) still using the old Dubna INC approximations 
that were obtained by Gudima et al. |18j 36 years ago, using the measurements available at 
that time. For instance, in the case of two-body reactions, the cosine of the angle of emission 
of secondary particles in the cm. system is calculated by the Dubna INC as a function of a 
random number ^, distributed uniformly in the interval [0,1] as 
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The coefficients a„fc were fitted to the then available experimental data at a number of incident 
kinetic energies Tj , then interpolated and extrapolated to other energies (see details in j^Q ^] 
and references therein). The distribution of secondary particles over the azimuthal angle (p is 
assumed isotropic. For elementary interactions with more than two particles in the final state, 
the Dubna INC uses the statistical model to simulate the angles and energies of products (see 
details in pUj). 

For the improved versions of our codes referred to as CEM03 and LAQGSM03, respectively, 
we use currently available experimental data and recently published systematics proposed by 
other authors to develop new approximations for angular and energy distributions of particles 
produced in nucleon-nucleon and photon-proton interactions. So, for pp, np, and nn interactions 
at energies up to 2 GeV, we did not have to develop our own approximations analogous to the 
ones described by Eqs. (1) and (2), since reliable systematics have been developed recently by 
Cugnon et al. for the Liege INC pHI) then improved further by Duarte for the BRIC code [20] : 
we simply incorporate into CEM03 and LAQGSM03 the systematics by Duarte [20]. Similarly, 
for interactions, we take advantage of the event generators for ^p and reactions from the 
Moscow INC 21, kindly sent us by Dr. Igor Pshenichnov. In our codes, we use part of a large 
data file with smooth approximations through presently available experimental data from the 
Moscow INC (21) and have developed a simple and fast algorithm to simulate unambiguously 
da I d£L and to choose the corresponding value of for any , using a single random number ^ 
uniformly distributed in the interval [0,1] Hj. For other elementary interactions, we fit new sets 
of parameters from Eq. (1) at different Tj for which we found data, then we approximated 
the energy dependences of the parameters a^k in Eq. (2) using the fitted values of a„. 

Examples of angular distributions of secondary particles from np and 7p reactions at several 
energies are shown in Figs. 1 and 2 of Ref. 0. The new approximations from CEM03 and 
LAQGSM03 reproduce the experimental data much better than the old Dubna INC used in our 
previous code versions (and in several other codes developed from the Dubna INC) and allow 
us to describe better particle spectra from different reactions on nuclei (see, e.g.. Fig. 3 in [5]). 

Third, we have improved the description of complex-particle spectra in CEM03 and 
LAQGSM03. This was done by refining the coalescence model used in our codes, by developing 
a better approach to estimate the probability of complex-particle emission at the preequilibrium 
stage of a reaction, and by incorporating into our codes the known systematics for angular 
distributions of complex particles developed by Kalbach (see details in jSj). 

Fig. 1 shows examples of proton, deuteron, and triton spectra from 542 MeV p + Cu and Bi 
calculated by the recently improved version of CEM compared with experimental data [221 ■ 
see that although the CEM03.01 results do not coincide exactly with these experimental data, 
the agreement is comparable to that provided by other modern models like FLUKA (see Fig. 
3 in |2n]) and the latest version of the Liege INC by Cugnon et al. (see Fig. 3 in [25), and is 
significantly better than that obtained a decade ago with the initial version of CEM JOl (see 
Figs. 6, 7, 9 and 10 in [2S1). 

We note that CEM03.01 also describes reasonably well complex-particle (and nucleon) spectra 
from reactions at energies below 100 MeV, where more sophisticated microscopic codes like 
GNASH j26j or TALYS [221 are usually used to produce data libraries for applications. For 
example, from Fig. 6 of Ref. [3] one can see that CEM03.01 describes the nucleon and complex- 
particles spectra from 62.9 MeV p + ^'^^Pb measured recently [2H| in Louvain-la-Neuve for the 
HINDAS project quite well and agrees with the data about as well as results by TALYS |27| . 
GNASH ,26 , FLUKA [2Sj, or INCL4 CTOl published in Figs. 20 and 21 of Ref. [^. 
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Figure 1. Measured |221 and calculated proton, 
deuteron, and triton spectra from 542 MeV p + 
Cu and Bi, respectively. 



Next we extended CEM03 to calculate 
reactions induced by both monochromatic 
and bremsstrahlung photons, as described 
in detail in [2l. As one can see from the 
example shown in Fig. 2 and from many 
other results published in ISHHl, CEM03.01 
now describes reasonably well photonuclear 
reactions at incident energies from about 30 
MeV to ~ 1.5 GeV. It can be used as well 
for estimation of photonuclear reactions for 
applications at higher energies, up to ~ 5 
GeV (see details in 1301). 
Finally, after making the improvements to 
the INC and preequilibrium parts GEM (and 
LAQGSM) as described above, the mean 
values of the mass and charge numbers, 
A and Z of the excited compound nuclei 
produced after the preequilibrium stage of 
nuclear reactions and their mean excitation 
energy E* have changed slightly, which 
affects the probability of heavy compound 
nuclei (especially preactinides) to fission. 
This means that the procedure of fitting 
the ratio of the level-density parameters of 
fissioning nuclei at the saddle point, af, and 
for evaporation from compound nuclei, a„. 




Product mass number A 



Figure 2. CEM03.01 results (green symbols) for 
the isotopic yields of products from Eq = 1 GeV 
bremsstrahlung-induced reactions on ^^^Au compared 
with experimental data (red symbols) |22j and calculations 
by PIGA3/GEM (blue symbols); the PICA3/GEM results 
are from several publications and are presented in Fig. 18 
of |29j . The black curves represent approximations based 
on experimental data [22] (see more details in 0). 
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Figure 3. Experimental [3Tj 
proton-induced fission cross sec- 
tions of I'^^W, ^'^^W, ^^^W, and 
182-^ compared with improved (red 
solid lines) and old (blue dashed 
lines, from |SI]) CEM03.01 calcu- 
lations. 



i.e., af/ttn which we performed in 30 to provide the best description by GEM2k and LAQGSM 

of fission cross sections is no longer correct. We redo for GEM03.01 and LAQGSM03.01 the 
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determination of the level-density ratio done in for CEM2k and LAQGSM, ensuring that 
the latest versions of our codes describe as well as possible fission cross sections from various 
reactions. Fig. 3 shows examples of fission cross sections for proton-induced reactions on ^^^W, 
^^^W, ^^^W and ^^^W. One can see that the improved CEM03.01 reproduces the recent Uppsala 
measurements jH^ of proton-induced fission cross sections. 

Our calculations have shown that CEM03.01 results also agree reasonably well with the recent 
Uppsala and Saint-Petersburg 35 data and older measurements at LANL by Parrish 

Staples et al. jSH] for neutron-induced fission cross sections. Results by LAQGSM03.01 for these 
and other fission cross sections practically coincide with ones by CEM03.01, as is expected from 
the fitting process [HUj . 

The initial version of LAQGSM [11, just like its precursor QGSM Tl", did not consider 
photonuclear reactions, while its 2003 version LAQGSM03 6 and CEM03.01 describe such 
reactions only for energies up to about 1.5 GeV. This is not convenient when using our codes 
to solve problems for PRAD, NASA, and other high-energy applications or to analyze future 
high-energy measurements at the Thomas Jefferson National Accelerator Facility (CEBAF), 
where photons with much higher energy are created and need to be simulated by an event 
generator in a transport code. To address this problem, we have extended LAQGSM03 to 
describe photonuclear reactions at energies up to tens of GeV. For this, we took advantage 
of the high-energy event generators for 'jp and jn elementary interactions from the Moscow 
high-energy photonuclear reaction model 21 kindly sent us by one of its co-authors, Dr. Igor 
Pshenichnov, as mentioned previously. We have incorporated into LAQGSM03 56 channels to 
consider 'jp elementary interactions during the cascade stage of reactions, and 56 channels for jn 
interactions. These reaction channels new to LAQGSM03 are listed in Table 1 of Ref. [01 together 
with all corresponding details. The improved version of our code extended in this manner is 
referred to as LAQGSM03.01; it allows us to describe photonuclear reactions at energies both 
below ~ 1.5 GeV where the earlier LAQGSM03 and CEM03.01 work (see, e.g. Fig. 11 in "W), 
and at higher photon energies. Figs. 4 and 5 show two examples of results by LAQGSM03.01 
for high-energy photonuclear reactions; more such results may be found in [21 . We note that 
to the best of our knowledge, we are able to describe with LAQGSM03.01 the data shown in 
Figs. 4 and 5 (and other similar reactions analyzed in ^) for the first time; we do not know of 
any publication or oral presentation where these measurements were reproduced by a theoretical 
model, event generator, or transport code. 

We have benchmarked the new versions of our codes, CEM03.01 and LAQGSM03.01, on a 
variety of particle-particle, particle- nucleus, and nucleus-nucleus reactions at energies from 10 
MeV to 800 GeV per nucleon and find that they describe reactions generally much better than 
their predecessors. Two examples of this work are shown in Figs. 6 and 7. We note that the 400 
GeV experimental data compared in Fig. 7 with our LAQGSM03.01 results and other similar 
measurements of K^, K~, p, d, t, ^He, and ^He spectra from ^^^Ta, as well as all measured 
spectra from Cu, Al, and C are described here simultaneously within a single approach 

for the first time: Though these data were measured at FNAL 25 years ago with a hope of 
revealing some "exotic" or unknown mechanisms of nuclear reactions leading to the production 
of the measured so called "cumulative" (i.e., kinematically forbidden for quasi- free intranuclear 
projectile-nucleon collisions) particles, we do not know any publication or oral presentation where 
these measurements were reproduced simultaneously within a single approach by a theoretical 
model, event generator, or transport code. It is noteworthy that LAQGSM03.01 describes quite 
well all cumulative particle spectra measured in in single approach, without any fitting 

or free parameters, and without involving any "exotic" reaction mechanisms. This is also true 
for the photonuclear cumulative proton yields shown in Fig. 4. These results do not imply, of 
course, that the proton- or 7-nucleus interaction physics is completely described by the reaction 
mechanisms considered by LAQGSM03.01. Our present results do not exclude some contribution 
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Figure 4. Proton spectra at 30, 60, 90, and 
150 degrees from interaction of bremsstrahlung 
7 quanta of maximum energy Eq = 4.5 GeV 
with 12c, 27 Al, ^^Cu, and ^ospb. Experimental 
values shown by symbols are from [SH EHl EHl 
while histograms show results by LAQGSM03.01. 
Similar results 9 were obtained by LAQGSM03.01 
at £^0 = 3.0 and 2.0 GeV, as well as for photopion 
spectra |^. To the best of our knowledge, 
using LAQGSM03.01 we are able to describe these 
data measured 28 years ago at Yerevan for the first 
time. 
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Figure 5. Detailed comparison 

between experimental yields jl2] and 
those calculated by LAQGSM03.01 of 
radioactive products from the interaction 
of bremsstrahlung 7 quanta of maximum 
energy 4.5 GeV with ^^Nb. The 
cumulative yields are labeled as "c" while 
the independent cross sections, as "i". 
To the best of our knowledge, using 
LAQGSM03.01, we are able to describe 
these data measured 24 years ago at 
Yerevan for the first time. 

to the production of these cumulative par- 
ticles from other reaction mechanisms not 
considered by LAQGSM03.01. But the 
contribution from "exotic" mechanisms to 
cumulative particle production from these 
high-energy reactions seems to be small; 
inclusive particle spectra are not sensitive 



enough for an unambiguous determination of the mechanisms of particle production, just as 
observed heretofore at intermediate and low energies "25" . 

Finally, we wish to note that CEM03.01 and LAQGSM03.01 and their predecessors have good 
predictive powers and describe various nuclear reactions similarly to and often better than other 
current models do; therefore they can be used as reliable event generators in transport codes for 
applications. Fig. 8 shows just one example to demonstrate this. This figure was made by adding 
the recent fission- fragment yield data from the reactions 1 GeV/A ^38^ _j_ (brown circles) 

to Fig. A2.3 of Ref. |Hj, which showed our calculations of several models. The calculations were 
done about two years before the data [^SI became available to us (the dates of all our calculations 
are shown in the legend of Fig. 8). The prediction by our LAQGSM+GEM2, a predecessor of 
LAQGSM03.01 for the fission-fragment mass distribution from this reaction agrees very well 
with the later measurements ^53, • The agreement of LAQGSM+GEM2 results with the older 
data ^52^ (magenta circles) on spallation product yields (which were available to us before the 
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Figure 6. Invariant tt~ yield per central 
Au + Au collision at 8 GeV/A as calculated 
by LAQGSM03.01 (histograms) and measured 
in [331 (symbols). Red solid histograms and 
open circles are forward production, while 
blue dashed histograms and open triangles are 
backward production. Midrapidity (the upper 
set) data is shown unsealed, while 0.1-unit 
rapidity slices are scaled down by successive 
factors of 10. 



calculations were performed) are also satis- 
factory and are comparable to that provided 
by LAHET3 HH calculations using: 

1) ISABEL I3H1 INC with the Dresner [IS] 
evaporation and RAL [SU] fission models; 

2) the Liege INC code INCL by Cugnon et 
al. merged with the Dresner [13] evapo- 
ration and RAL (50j fission models; 

3) INCL HH] merged with ABLA/PROFI 
fission/evaporation model by Schmidt et al. 

m- 

4. Applications 

CEM03.01 and LAQGSM03.01 and their pre- 
decessors are used as stand-alone codes to 
study different nuclear reactions for applica- 
tions and fundamental nuclear physics (see, 
e.g., [31 [3 [HI [ISl and references therein). 
We outline below only two examples of ap- 
plications. 

The first application of CEM97 and 
CEM95, predecessors of CEM03.01, we like 
to mention here is for a medical isotope pro- 
duction study. Several years ago, a detailed 
study of the production of 22 isotopes for 
medical and industrial applications by high- 
energy protons and neutrons at the Acceler- 
ator Production of Tritium Facility project 
(APT) was performed (see IMl and ref- 
erences therein). The production rate of a 
radioisotope can be obtained from the inte- 
gral of the flux and cross section leading to 
the direct production of the radioisotope as 



a reaction product. Additional production is realized from other radionuclides that decay 
to the desired product. Evaluation of production rates requires knowledge of the neutron 
and proton fluxes at some position in the production facility and cross sections leading to 
production of the desired radionuclide and its progenitors. We used MCNPX version 2.1.1 
to calculate neutron and proton fluxes throughout the APT model, while cross sections for 
reactions most likely to lead to our desired products were evaluated using all available to us 
experimental data and calculations by CEM95 and CEM97 at energies above 100 MeV, and the 
"150 MeV" activation library calculated by M. B. Chadwick with the HMS-ALICE code [S3 
supplemented by the European Activation File EAF-97, Rev. 1 with some recent improvements 
by M. Herman (see ^j] and references therein), at lower energies. Our 684 page detailed 
report on this study '55 , with 37 tables and 264 color figures is available on the Web under: 
: / /t2. lanl.gov /publications /publications. html 
A second application of our codes is for astrophysics. A successful model for cosmic-ray 
propagation in the Galaxy used nowadays to study different challenging questions of astrophysics 
is realized in the code GALPROP by Strong and Moskalenko (see [SHI and references therein). 




p (GeV/c) p (GeV/c) 

Figure 7. Invariant spectra of p, d, vr"'", 
and vr^ from tiie reaction 400 GeV p + ^'^^Ta. 
Experimental data for p are from [llj, for d from 
|45j ■ and for pions from Similar results 

are obtained by LAQGSM03.01 for the measured 
, K~ , p, d, t, ^He, and ^He spectra from 
^^^Ta, as well as for all measured spectra from 
Cu, Al, and C [m ESI EE]. To the best of our 
knowledge, using LAQGSM03.01 we are able to 
describe simultaneously within a single approach 
all these data measured 26 years ago at FNAL 
for the first time. It is interesting that we do 
not need to consider any "exotic" nuclear-reaction 
mechanisms to describe reasonably well all these 
measured "cumulative" particle spectra. 



So, GALPROP was used recently to 
study propagation of different particles 
and nuclei in the Galaxy and to improve 
estimates of the Galactic halo size (see, 
e.g., 0, [nni-IMl and references therein). 
However, GALPROP relied initially on 
nuclear reaction cross sections calculated 
using the phenomenological systematics 
by Webber et al. [HSl and Silberberg et al. 
|64j and on scarce available experimental 
data, as is usually done in astrophysics. 
As an example. Fig. 9 shows only one 
excitation function for the production 
of ^^Al from the reaction p + "■'^*Si 
used by GALPROP along with hundreds 
of other excitation functions. We see 
that the phenomenological systematics 
inSl inn do not reproduce correctly the 
available experimental cross sections for 
this reaction. As one can see from 
Fig. 10, using in GALPROP cross 
sections provided by phenomenological 
systematics [HSl IHl] leads to quite big 
uncertainties in the derived Galactic halo 
size (3-7 kpc), while using evaluated cross 
sections based on available experimental 
data and calculations by CEM2k 
reduces uncertainties and limits the 
derived Galactic halo size to only 4-6 
kpc. This is why we have used several 
versions of our GEM and LAQGSM 
codes to calculate cross sections of 
interest to astrophysics, used thereafter 
together with available experimental data 
to produce a number of reliable evaluated 
excitation functions for astrophysical 
needs (see, e.g., jHS] and 

references therein). 



Finally, we mention that different versions of our GEM and LAQGSM codes are incorporated 
wholly as event-generators, or in part in different transport codes used in applications, 
among them CASCADE MARS 67, MCNPX 68 GEANT4 m, SHIELD 
RTSfcT jZH, SONET [Z2], CALOR 73 , HETC-3STEP [ZH, CASCADE/INPE 
HADRON 76 , CASCADO [ZZj, CAMO 78^ and others. Thus our codes are automatically 
employed in various nuclear applications where these transport codes are used. 

5. Summary 

Improved versions of the cascade-exciton model of nuclear reactions and of the Los Alamos 
quark-gluon string model have been developed recently at LANL and implemented in the codes 
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CEM03.01 and LAQGSM03.01 as event-generators for transport codes MCNP6 (2^!, MCNPX 
|68j . and MARS 67 . Our codes were previously incorporated into MARS and are now being 
incorporated into MCNPX and MCNP6. CEM03.01 was made available to the public via RSICC 
at Oak Ridge. We also plan to make LAQGSM03.01 available to the public via RSICC in the 
future. 
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Figure 8. Comparison of measured mass (left panel) and charge (right panel) distributions of 
the nuclides produced in the reactions 1 GeV/A 238 ^ _j_ ^ -^^ith results from LAHET3 ^47, using 
the IS ABEL+Dresner/ Atchison |lHl-|Sni (green lines), INCL+Dresner /Atchison ISU] 
(blue line), and the INCL+ABLA [H EI] (red lines) models, and by our LAQGSM+GEM2 
(maroon lines), respectively. Experimental data for spallation products shown by magenta 
circles are from |52| and were available to us prior to the calculations being done at the dates 
indicated in the legend, while the fission-product data [SSI; shown here by brown circles, became 
available to us only one year after all the calculations were published in Fig. A2.3 of Ref. [S]. 
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